Phys. Rev. C (1997) in press. 



Proton and Antiproton Distributions at Mid-Rapidity in 
Proton-Nucleus and Sulphur-Nucleus Collisions 

I.G. Bearden 1 , H. B0ggild 1 , J. Boissevain 2 , J. Dodd 3 , S. Esumi 4 ' a , C.W. Fabjan 5 , 
D.E. Fields 2 ' 6 , A. Franz 5 > c , A.G. Hansen 1 , E.B. Holzer 5 , T.J. Humanic 6 , B.V. Jacak 2 - d , 
R. Jayanti 6 ' 7 , H. Kalechofsky 7 , Y.Y. Lee 7 , M. Leltchouk 3 , B. Lorstad 8 , N. Maeda 4 ' 6 , 
A. Medvedev 3 , A. Miyabayashi 8 , M. Murray 9 , S. Nishimura 4 ^, S.U. Pandey 7 ' 9 , F. Piuz 5 , 
V. Polychronakos 10 , M. Potekhin 3 , G. Poulard 5 , A. Sakaguchi 4 '' 1 , M. Sarabura 2 , 
M. Spegel 5 , J. Simon-Gillo 2 , W. Sondheim 2 , T. Sugitate 4 , J.P Sullivan 2 , Y. Sumi 4 , 
H. van Hecke 2 , W.J. Willis 3 , K. Wolf 9 '*, and N. Xu 2j 

(The NA44 Collaboration) 
1 Niels Bohr Institute, DK-2100 Copenhagen, Denmark. 
2 Los Alamos National Laboratory, Los Alamos, NM 87545, USA. 
3 Columbia University, New York, NY 10027, USA. 
4 Hiroshima University, Higashi- Hiroshima 724, Japan. 
5 CERN, CH-1211 Geneva 23, Switzerland. 
6 Ohio State University, Columbus, OH 43210, USA. 
7 University of Pittsburgh, Pittsburgh, PA 15260, USA. 
8 University of Lund, S-22362 Lund, Sweden. 

9 Texas A&M University, College Station, TX 77843, USA. 

10 Brookhaven National Laboratory, Upton, NY 11973, USA. 

a Now at Heidelberg University, Heidelberg, D-69120, Germany 
b Now at University of New Mexico, Albuquerque, NM 87185 
c Now at Brookhaven National Laboratory, Upton, NY 11973 
d Now at State University of New York Stony Brook, Stony Brook, NY 11794 
e Now at Florida State University, Tallahassee, FL 32306 

1 



f Now at University of Tsukuba, Ibaraki 305, Japan 
g Now at Wayne State University, Detroit, MI 48201 
h Now at Osaka University, Osaka 560, Japan 
1 deceased, 

i Now at Lawrence Berkeley National Laboratory, Berkeley, CA 94720 



Abstract 



Experiment NA44 has measured proton and antiproton distributions at 
mid-rapidity in sulphur and proton collisions with nuclear targets at 200 and 
450 GeV/c per nucleon respectively. The inverse slopes of transverse mass 
distributions increase with system size for both protons and antiprotons but 
are slightly lower for antiprotons. This could happen if antiprotons are an- 
nihilated in the nuclear medium. The antiproton yield increases with system 
size and centrality and is largest at mid-rapidity. The proton yield also in- 
creases with system size and centrality, but decreases from backward rapidity 
to mid-rapidity. The stopping of protons at these energies lies between the 
full stopping and nuclear transparency scenarios. The data are in reason- 
able agreement with RQMD predictions except for the antiproton yields from 

sulphur- nucleus collisions. 

PACS numbers: 25.75.-q 13.85.-t 13.60.Rj 
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I. INTRODUCTION 



Nucleus-nucleus collisions at ultrarelativistic energies create hadronic matter at high 
energy density. The distributions of baryons at mid-rapidity provide a sensitive probe of the 
collision dynamics. In particular, the stopping power determines how much of the incoming 
longitudinal energy is available for excitation of the system. These collisions have been 
described by microscopic models incorporating hadron production and rescattering, see for 
example [IJ. The data reported here impose constraints on the amount of stopping of baryons 
in such models. 

Enhanced production of antibaryons may indicate formation of a state of matter in which 
the quarks and gluons are deconfined . Such enhancement may be hidden, however, 

by antibaryon annihilation with baryons ; the antiproton survival probability is sensitive 
to both the collision environment and the antiproton formation time. The antiproton and 
proton distributions may also reflect the degree of thermalization achieved and, by comparing 
distributions from light and heavy systems, allow detailed studies of rescattering. 

We present proton and antiproton measurements using the NA44 spectrometer from 
pBe (to approximate pp), proton-nucleus and nucleus-nucleus interactions. This allows a 
systematic study as a function of the size of the central region and different conditions in 
the surrounding hadronic matter. 

These systematic studies are aided by use of an event generator. The RQMD model, 
version 1.08 [1]], is a microscopic phase space approach, based on resonance and string exci- 
tation and fragmentation with subsequent hadronic collisions. RQMD includes annihilation 
of antiprotons in the hadronic medium when they collide with baryons J7|. We study a 
feature of RQMD called color ropes by the authors of the model. RQMD uses a string 
model of particle production from each nucleon-nucleon collision. In a heavy-ion collision, 
where there are numerous nucleon-nucleon collisions, the density of these strings is high 
and multiple strings overlap. Overlapping strings do not fragment independently but form 
'ropes', chromoelectric flux-tubes whose sources are color octet charge states rather than 
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the color singlet charges of normal strings ||. These ropes represent a collective effect in 
nucleus-nucleus collisions, and have been shown to enhance both strangeness and baryon 
pair production ||. 

II. EXPERIMENT 

The NA44 experiment is shown in Figure [l|. Three conventional dipole magnets (Dl, D2, 
and D3) and three superconducting quadrupoles (Ql, Q2, and Q3) analyze the momentum 
and create a magnified image of the target in the spectrometer. The magnets focus particles 
from the target onto the first hodoscope (HI) such that the horizontal position along the 
hodoscope gives the total momentum. Two other hodoscopes (H2, H3) measure the angle 
of the track. The momentum acceptance is ± 20% of the nominal momentum setting. The 
angular coverage is approximately -5 to +78 mrad with respect to the beam in the horizontal 
plane and ± 5 mrad vertically. Only particles of a fixed charge sign are detected in a given 
spectrometer setting. Four settings are used to cover the mid-rapidity region in the pt range 
to 1.6 GeV. Figure shows the acceptance of the spectrometer in the y — pr plane for the 
4 and 8 GeV/c momentum settings when the spectrometer axis is at 44 and 131 mrad with 
respect to the beam. 

For the sulphur-nucleus data, the beam rate and time-of-flight start are determined using 
a Cherenkov beam counter (CX), with time resolution of approximately 35 ps f9|. For the 
proton-nucleus data, a scintillator counter is used to measure the beam rate. A second 
scintillator (To) is used to trigger on central events in sulphur-nucleus collisions by requiring 
a large pulse height (high charged particle multiplicity). The pseudorapidity coverage of T 
is roughly 1.3 to 3.5. For proton beams, To provides the interaction trigger by requiring 
that at least one charged particle hit the scintillator, and also provides the time-of-flight 
start with a time resolution of approximately 100 ps. A silicon pad detector measures 
the charged-particle multiplicity with 27r azimuthal acceptance in the pseudorapidity range 
1.8 < T] < 3.3. 
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The three scintillator hodoscopes (HI, H2 and H3) are used to track the particles and 
are divided into 50, 60 and 50 slats, respectively. The hodoscopes also provide time-of-flight 
with a time resolution of approximately 100 ps; particle identification relies primarily upon 
the third hodoscope. Two Cherenkov counters differentiate kaons and protons (CI: freon-12 
at 1.4 or 2.7 atm, depending on the spectrometer setting), and reject electrons and pions 
(C2: nitrogen/neon mixture at 1.0 or 1.3 atm). An appropriate combination of CI and 
C2 is used for each spectrometer setting to trigger on events with no pions or electrons in 
the spectrometer. Particles are identified by their time-of-flight, in combination with the 
Cherenkov information. Figure |3| illustrates the particle identification after pions have been 
vetoed by the Cherenkovs; kaons and protons are clearly separated. More details about the 



spectrometer are available in [10 



III. DATA ANALYSIS 

The proton and antiproton data samples after particle identification and quality cuts are 
shown in Table |. Also shown for each data set is the target thickness and the centrality, 
expressed as a fraction of the total inelastic cross-section. 

Tracks are reconstructed from the hit positions on the three hodoscopes, constrained 
by straight-line trajectories after the magnets. In order to construct the invariant cross- 
section, the raw distributions are corrected using a Monte Carlo simulation of the detector 
response. Simulated tracks are passed through the full analysis software chain and used 
to correct the data for geometrical acceptance, reconstruction efficiency and momentum 
resolution. Particles are generated according to an exponential distribution in transverse 
mass, ttit = \Jpt + m \i with the coefficient of the exponent determined iteratively from the 
data. 

The Cherenkovs reject pions with an efficiency of greater than 98% at the trigger level. 
Further offline selection reduces the pions to a few percent of the kaons. After time-of-flight 
and Cherenkov cuts, the residual kaon contamination of the proton sample is less than 3%. 
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The invariant cross-sections are presented as a function of ttlt - m p , where m p is the mass 
of the proton. The absolute normalization of each spectrum is calculated using the number of 
beam particles, the target thickness, the fraction of interactions satisfying the trigger, and the 
measured live-time of the data-acquisition system. For the SA data, the centrality selection 
is determined by comparing the pulse height distribution in the T counter for central and 
unbiased beam triggers. For pA systems, the fraction of inelastic collisions producing at 
least one hit in the interaction (To) counter is modeled with the event generators RQMD 
U and FRITIOF []ll , E| . The errors on the centralities for the pA data in Table H reflect 



the systematic uncertainty on this fraction from comparing the two models. The pA data 
are also corrected for the efficiency of the T counter. The resulting centrality fractions are 
indicated in Table |. 

The proton-nucleus data are corrected for non-target background. The largest corrections 
are 7.4% and 6.7% on the absolute cross-sections for protons and antiprotons from pBe 
collisions. This correction does not affect the shape of the distribution. The corrections to 
the nucleus-nucleus data are negligible. The cross-sections are also corrected for the proton 
identification efficiency, and for the effects of selecting events with no accompanying pion 
or electron. The pion (electron) veto correction is determined from the number of protons 
in runs for which the pions (electrons) are not vetoed. The effect of hadronic interactions 
of the produced particles in the material of the spectrometer has been studied using a 
detailed GEANT simulation. These interactions, including annihilation of antiprotons in the 
spectrometer material, do not distort the shape of the measured transverse mass distributions 
but result in a reduction in the observed yields of about 11% for protons and 17% for 
antiprotons. The data are corrected for these losses. 

The invariant cross-sections, measured in the NA44 acceptance, are generally well de- 
scribed by exponentials in transverse mass (see Equation [1]). The proton and antiproton 
rapidity densities (dN/dy) are calculated by integration of the normalized my distributions, 
with the fitted coefficient of the exponent (the 'inverse slope') used to extrapolate to high 
rriT, beyond the region of measurement. The statistical error on this extrapolation is cal- 

6 



culated using the full error matrix from the fit of Equation [1] (see Section |VI|) to the itlt 
spectrum. The corresponding systematic error is included in Table y. 

IV. SYSTEMATIC ERRORS 

The momentum scale of the spectrometer is verified with a second, independent measure- 
ment of the momentum using the multiwire proportional chambers (MWPCs 1-4 in Figure 
[l]) and dipole magnet (D3), yielding a systematic error of 1.6% on the px scale. Since the 
NA44 spectrometer has acceptance for both positive and negative px, a systematic offset 
in px can be checked by requiring symmetry around px = 0. The resulting uncertainty 
on the origin of the pt scale is 7 MeV/c for the 8 GeV setting. A detailed Monte Carlo 
simulation, including multiple scattering and detector granularity, is used to correct for the 
finite resolution of the spectrometer, and introduces a systematic error in p? of 0.15%. 

Systematic errors on the inverse slopes of the transverse mass distributions are estimated 
by comparing the inverse slope determined from the 131 mrad data to the inverse slope de- 
termined from both angle settings, and are less than 5%(15%) for the 8 GeV(4 GeV) setting 
data. The systematic errors due to the spectrometer acceptance correction are estimated 
from the sensitivity of the extracted slope to the fit range used, and by measuring the slope 
determined from the ratio of the cross-sections corresponding to the 'central ray' of the 
spectrometer at both the 131 and 44 mrad settings. In this ratio the acceptance corrections 
cancel since the particles have the same path through the spectrometer. The total uncer- 
tainty is 10 MeV/c for the mid-rapidity (y— 2.3-2.9) data, and 10-20 MeV/c for the lower 
rapidity (y=1.9-2.3) data. 

The error in the absolute normalizations is dominated by the uncertainty in the fraction 
of the total cross-section selected by the NA44 centrality trigger, and by the pion and 
electron veto corrections. The relative error in the centrality is 6% for both the SS and SPb 
data, resulting in a systematic uncertainty of 6% in dN/dy. For the proton-nucleus data, 
the trigger bias is determined by modeling the acceptance of the T counter, as described 



7 



above. The resulting dN/dy values are sensitive to the charged particle distribution from 
the two models, giving an uncertainty of 1.5% for the pBe data, increasing to 3% for pPb. 
Corrections for the fraction of events vetoed by pions are significant only for the sulphur- 
nucleus data at the 44 mrad setting. The uncertainties in these corrections are 10% for 
SPb collisions and 5% for SS collisions. Additional errors on the dN/dy values arise from 
a 5% uncertainty in the determination of the data-acquisition live-time, and a 5% beam 
rate dependence of the pion veto correction for the SPb data. For the pA data, there is 
an additional 6% error due to the correction for the efficiency of the interaction counter. 
The correction for non-target background contributes a negligible systematic error to the 
absolute cross- sect ions for proton-nucleus data. The total systematic errors on the measured 
inverse slopes and dN/dy values are given in Table ||. 



V. FEED-DOWN FROM WEAK DECAYS 

Weak decays of strange baryons are a significant source of protons and antiprotons, 
and contribute to the yields measured in the NA44 spectrometer. A strange baryon may 
travel several centimeters from the target before decaying weakly to a proton and a pion. 
Such a proton may be reconstructed in the spectrometer with the wrong momentum. The 
sensitivity of the data to this feed-down has been studied using a GEANT simulation of 
the spectrometer with particle distributions and yields taken from RQMD. Of the strange 
baryons, only A and E + decay weakly to protons (with the corresponding antiparticles 
decaying to antiprotons). Heavier strange baryons also contribute via sequential decays. 

The fraction of the measured protons (antiprotons) arising from these decays is calculated 
as a function of rapidity and transverse momentum within the NA44 acceptance, and used to 
determine 'feed-down factors' for the rriT and dN/dy distributions. These are multiplicative 
factors which provide an estimate of the contribution of A and E decays to the measured 
distributions, and could be applied to the numbers in Tables [TV] and |V| to estimate the inverse 
slopes and yields of 'direct' protons (antiprotons). Figure [| shows the effect of feed-down 
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from weak decays on the RQMD proton and antiproton transverse mass distributions for 
SPb collisions. Including A and E decays tends to make the distributions steeper since the 
protons arising from weak decays contribute more at low p?. Table [TII| shows the feed-down 
factors for the inverse slopes and yields of the data due to A and E decays. The values 
are calculated using the A/p and E/p ratios from RQMD. The 'errors' reflect the result of 
increasing and decreasing the (A + E)/p ratio in the model by a factor of 1.5. This factor of 
1.5 is consistent with the scale of the discrepancies in the published data on A production in 
nucleus- nucleus collisions. For central SS collisions, RQMD is consistent with the A/p ratio 
measured by NA35 [O]. However for SPb, there is a factor of 2 discrepancy between the 



measurement of A production by NA36 []14j and the scaled SAg data of NA35. The RQMD 
prediction lies between the results of the two experiments. The feed-down factors have 
been calculated explicitly for protons from pBe, pPb and SPb using the complete GEANT 
simulation, and are scaled according to the respective RQMD (A + E)/p ratios for protons 
from pS and SS, and for antiprotons from all systems. A factor of 1.5 variation in the 
(A + E)/p ratio is also assumed in calculating the antiproton factors. The dependence 
of the feed-down factors is mainly determined by the experimental acceptance and not by 
the tut dependence of the (A + S)/p ratio from RQMD. 

As these factors necessarily contain some dependence on models they have not been 



applied to the data but are listed in Table 0TI] so that the reader can appreciate the sensitivity 
of the data to A and E decays. 

VI. RESULTS 

The invariant cross-sections for protons and antiprotons from pBe, pS, pPb, SS and 
SPb collisions are shown in Figures and ^ as a function of - m p . The transverse 
mass distributions are generally well described by exponentials in the region of the NA44 
acceptance: 

1 Ed3(7 = Ce -{m T -m p )/T /j\ 

cr dp 3 



where C is a constant and T the inverse logarithmic slope. The inverse slope parameters 
obtained by fitting the proton and antiproton data to Equation [I] are given in Table [TV], 
and plotted in Figure |7|. The inverse slopes for both protons and antiprotons increase with 
system size. For protons, the inverse slopes are higher at mid-rapidity (y— 2.3-2.9) than at 
more backward rapidities (y =1.9-2. 3). This effect is not seen for antiprotons, where the 
inverse slopes are similar in the two rapidity intervals. Note that the errors on the backward 
rapidity data are significant. The inverse slopes for antiprotons are generally somewhat lower 
than for protons at mid-rapidity, but are comparable in the backward rapidity interval. 

The proton and antiproton rapidity densities (dN/dy) are listed in Table [V| and plot- 
ted in Figure |8|. Proton yields increase with system size, and are significantly larger in 
nucleus-nucleus collisions than in proton-nucleus collisions. More protons are produced in 
the backward rapidity interval (y =1.9-2. 3) than at mid-rapidity (y=2.3-2.9), particularly for 
proton-nucleus collisions. The antiproton yields are lower than the proton yields and grow 
less rapidly with increasing system size: the increase in the antiproton yield between pBe 
and pPb is less than 50%. Comparing antiproton yields in the two rapidity intervals, there is 
essentially no difference for proton-nucleus collisions. In nucleus-nucleus collisions however, 
antiproton production is notably smaller backwards of mid-rapidity. 

Figure || shows the ratio of p to p yields for the various projectile-target systems. Note 
that the systematic errors described in Table || cancel in this ratio. The p/p ratio decreases 
by a factor of 4 from pBe to SPb. Most of this decrease with system size occurs between 
pPb and SS. The target dependence of the ratio is stronger in sulphur-nucleus than in 
proton-nucleus collisions. Comparing the two rapidity intervals, the p/p ratio is larger at 
mid-rapidity in all cases. 

The beam momentum for the pA data is 450 GeV/c, corresponding to a beam rapidity 
of 7, while for the SA data the beam momentum is 200 GeV/c per nucleon, corresponding 
to a beam rapidity of 6. The energy dependence of proton and antiproton production in 
pp collisions has been studied at the ISR |T5[. Using the parametrizations of the proton 



and antiproton cross-sections as a function of center-of-mass energy from ||15|| , the effect of 

10 



the different beam momenta for the pA and AA data on the systematic behaviour of the 
p/p ratio can be estimated. Decreasing the beam momentum from 450 to 200 GeV/c per 
nucleon, this energy scaling implies that the p yield decreases by 12 ± 1% and the p yield 
decreases by 1.5 ± 0.1%. Thus the decrease in beam energy cannot explain the decrease of 
the p/p ratio. Rather it reflects the fact that at mid-rapidity most protons are not produced 
but originate from the target or projectile. 

In order to study the centrality dependence of proton and antiproton production, the SPb 
data are divided into two subsamples containing the 11-6% and 6-0% most central collisions 
respectively. The inverse slopes show no centrality dependence between the two subsamples. 
Figure [10] shows the proton and antiproton yields for these two different centrality selections. 
Production of both protons and antiprotons increases with centrality, although the yield of 
protons rises faster. 



VII. DISCUSSION 

The precision of the data and the range of the systems studied provide strong constraints 
on models of proton and antiproton production, rescattering and annihilation. We have 
compared our data in detail to the RQMD model. Figure [lT] compares the transverse mass 
distributions of protons and antiprotons from SPb collisions to predictions from RQMD with 
and without rope formation. The contribution of feed-down from A and £ decays is included 
in the model predictions, determined using the GEANT simulation described above. The 
shapes of the spectra are generally reproduced by the model. This is true also for the lighter 
collision systems. The RQMD distributions are then used to determine inverse slopes and 
yields from the model in the same way as for the data. 

Comparing different collision systems, the inverse slopes of both protons and antipro- 
tons increase with system size. The relative increase, from pBe to SPb, is similar in both 
cases. NA44 has previously reported an increase in the inverse slopes of kaons, protons 



and antiprotons produced at mid-rapidity in symmetric collision systems JT8I]. These data 
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extend this trend to asymmetric systems, and to more backward rapidities (|/=1.9-2.3). For 
symmetric systems, the inverse slopes of protons and antiprotons are equal, while for the 
asymmetric systems pS,pPb and SPb at y=2.3-2.9 the inverse slopes of protons are higher 
than those of antiprotons. 

The increase of the proton inverse slope with system size, which continues through PbPb 



collisions ||18|| , is a result of the increasing number of produced particles and consequent 
rescattering. The large number of secondary collisions causes the hadronic system to ex- 
pand and the velocity boost from this collective expansion is visible as an in- 



creased inverse slope in the proton ttit distributions. One would expect this effect to be 
concentrated at mid-rapidity, and this is supported by the data which show that the inverse 
slopes of protons are higher at y=2.3-2.9 than at y=1.9-2.3. Though the antiproton inverse 
slope is also increased by secondary collisions, the RQMD model predicts that the observed 
antiprotons have suffered fewer rescatterings than the observed protons. This is because sec- 
ondary collisions with baryons can annihilate the antiprotons, and the number of baryons 



at mid-rapidity in the nucleus-nucleus collisions is substantial [IB]. This may explain why 
the inverse slopes of antiprotons do not increase from ?/=1.9-2.3 to y=2.3-2.9. 

RQMD reproduces the trend of increasing inverse slope with system size (Table [TV] and 
Figure |7|), and is in reasonable agreement on the absolute values of the inverse slopes, with 
the exception of the backward rapidity (y=l. 9-2.3) protons from sulphur-nucleus collisions, 
where the model predicts larger inverse slopes than are measured experimentally. For SS 
and SPb at y =1.9-2. 3, the RQMD proton distributions tend to deviate from exponentials 
in transverse mass, showing larger inverse slopes at low my (within the NA44 acceptance) 
than at high tut- For RQMD, there is no significant difference in the inverse slopes if rope 
formation is turned off. This might seem surprising for antiprotons since if rope formation 
is included, almost all antiprotons originate from ropes, and the production mechanisms of 
antiprotons from ropes and from conventional sources are very different. However, rescatter- 
ing would tend to mask any difference in the initial ttit distribution of antiprotons produced 
from either ropes or conventional sources. 
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The yields of protons and antiprotons both increase with system size, the proton yield 
rising faster. The increase of both yields and inverse slopes with target size is much stronger 
in SA than in pA collisions since target nucleons may be struck by more than one projectile 
nucleon. The p/p ratio decreases with system size, and is lower towards target rapidity than 
at mid-rapidity. This implies that most of the protons at mid-rapidity are not produced in 
the collision, but are remnants of the initial nuclei. The fraction and absolute numbers of 
such 'original' protons decreases from y=1.9-2.3 to y=2.3-2.9, indicating that the stopping 
of protons at these energies is incomplete. For pA collisions, the antiproton yield is essen- 
tially constant from y=1.9-2.3 to y=2.3-2.9, whereas for SA collisions more antiprotons are 
produced at y=2.3-2.9. 

For the SA data, no change in the inverse slopes of protons and antiprotons with cen- 
trality is seen within the 11% most central SPb, and 8% most central SS, collisions. Both 
proton and antiproton yields increase with centrality. The increase is larger for protons than 
for antiprotons, and both particles show a larger increase at y=2.3-2.9 than at y=1.9-2.3. 

RQMD predictions for proton and antiproton yields, with and without color ropes, are 
shown in Table [VT| and Figure [S[ The proton yields are generally well described by the model 
in both pA and SA collisions, although there is a tendency for RQMD to overpredict the 
yield in the more backward rapidity interval for sulphur-nucleus collisions. Antiproton yields 
from pA collisions are also in reasonable agreement with RQMD, however the data show a 
somewhat weaker target dependence than the model. The antiproton yield in SS collisions 
is consistent with the RQMD prediction without ropes, whereas the SPb yield is between 
the rope and no-rope predictions. With no rope formation, it is possible that RQMD could 
reproduce the observed SA antiproton yields with a smaller annihilation cross-section in the 
model. If, however, the rope hypothesis is valid, then either RQMD overpredicts antiproton 
production in nucleus-nucleus collisions, or it underpredicts the subsequent annihilation of 
the produced antiprotons in the surrounding medium. 

Preliminary results from NA44 on the 7r + yields near mid-rapidity (y=3. 0-4.0) show 
that the pion yield increases by a factor of 36 ± 4 from pBe to SPb while the proton 
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and antiproton yields increase by factors of 75 ± 4 and 20 ± 2 respectively. The fact 
that antiproton production increases less rapidly with system size than pion production 
naively suggests that the antiproton yield from SPb collisions may be lowered because of 
annihilation. 

NA35 has published data on "net protons" (i.e. the difference between protons and 
antiprotons) and antiprotons forpS", pAu, SS (3% most central) and SAu (6% most central) 
collisions [[i6i[T7| . In those cases for which comparison data are available, the inverse slopes 



measured by the two experiments are consistent. Figure [T^ shows a comparison of the 
rapidity densities for net protons and antiprotons from NA44 and NA35. There is good 
agreement between the two experiments. 



CONCLUSIONS 

These data constitute the first systematic measurement of proton and antiproton yields 
and spectra from pBe to SPb. The stopping of protons is incomplete at 200GeV/A. For 
SA collisions, target nucleons tend to be struck by more than one projectile nucleon: this 
causes the target dependence of the inverse slopes and yields to be stronger in SA than in 
pA collisions. As the size of the system increases, the increasing density of particles at mid- 
rapidity causes the protons to recscatter more often and so their mean tut, or inverse slope, 
increases. This mechanism is less efficient for antiprotons since they may annihilate when 
they rescatter. In spite of annihilation, the yield of antiprotons increases with system size 
and centrality: this increase is strongest at mid-rapidity. Comparisons with RQMD imply 
that antiproton annihilation is overestimated in the model or that some new mechanism is 
needed to account for antiproton production in sulphur-nucleus collisions. 
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TABLES 



System 


Centrality 


Angle 


Target 


y— 


1.9-2.3 


y=2.3-2.9 






(mrad) 


Thickness 


p 


P 


p 


p 


pBe 


84±2% 


44 


3.3% 


22594 


14925 


35630 


17880 






131 


3.3% 


11200 


3269 


5568 


1823 


pS 


90±2% 


44 


3.3% 


44505 


12081 


1609 


5392 






131 


3.3% 


51341 


1956 






pPb 


97±3% 


44 


4.7% 


14094 


13754 


1913 


4467 






131 


9.9% 


22908 


856 


13907 


2274 


SS 


8.7±0.5% 


44 


6.6% 


11135 


3200 


5960 


16660 






131 


6.6% 


18652 


2703 


18379 


3762 


SPb 


10.7±0.6% 


44 


5.9% 


13833 


2081 


12062 


7293 






131 


5.9% 


38178 


4528 


43261 


7299 



TABLE I. Centrality, target thickness and number of events for each spectrometer setting. The 
target thickness is quoted in nuclear collision lengths for the given system. 
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^ VCt (itTl 
kj Volclll 


Error on Inverse Slope 


HjL I Ul Oil (_l±^l / t.J. V 




y=1.9-2.3 


y=2.3-2.9 




pBc 


10 MeV/c 


10 MeV/c 


y /o 


pS 


10 MeV/c 


10 MeV/c 


9% 


P Pb 


10 MeV/c 


10 MeV/c 


10% 


SS 


20 MeV/c 


10 MeV/c 


9% 


SPb 


20 MeV/c 


10 MeV/c 


14% 



TABLE II. Systematic errors on the inverse slopes and dN/dy. 
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System 


Parameter 


y 

Proton 


■i r\ o o 

= 1.9-2. «i 

Antiproton 


y 

Proton 


o o o 

— 2.6-2. \3 

Antiproton 


pBe 


1 


i ni i n m 

1.(J1±U.U1 




1 no i n ni 

l.Uz±U.l)l 


1 no i n m 


ai\ / ay 


n nn-r-n hq 
U.yU^U.Uo 




n nn-r-n no 

u.yu^u.uo 


n Qc^-r-n r\A 
11.00^11.114 


pS 


l 


1 m _i_n m 
1.U1±U.U1 


1 no_i_n m 
i.Uz±U.Ui 


1 no_i_n no 
l.U3±U.l)z 


1.U7 ±U.U3 


aiN / ay 




n ^n-r-n nr 


n m -i-n r\A 


n QQ-r-n n^i 
U.oo^U.UO 


pFb 


1 


1 nn i n no 

i.uy±u.u2 


1 on i n o 

i.zy±u.z4 


i n a i n n i 

l.U4±U.Ul 


1 1 i i n no 

l.ll±U.Uo 


aiN / ay 




U. 1 z^pU.Ub 


n nQ-r-n no 
U.yo^U.Uz 


n qi -i-n n/i 
U. 01^11.114 


ss 


T 


1.05±0.01 


1.09±0.03 


1.02±0.01 


1.08±0.02 


dN/dy 


0.78=f0.05 


0.74=f0.06 


0.84=f0.04 


0.71=f0.06 


SPb 


T 


1.01±0.01 


1.09±0.02 


1.10±0.02 


1.12±0.03 


dN/dy 


0.77T0.06 


0.67T0.07 


0.82^0.05 


0.73T0.06 



TABLE III. 'Feed-down factors' for A and £ decays. T is the inverse slope of the itit spectrum. 
These are multiplicative factors which could be applied to the numbers in Tables [V and [V] to 
estimate the inverse slopes and yields of 'direct' protons (antiprotons). The errors reflect the result 
of increasing and decreasing the (A + ratio in RQMD by a factor of 1.5. The errors on T and 
dN/dy are anticorrelated. 
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y 


Fit Range 


System 


P 




P 


Ivl^lVllJ p 




^Vjt3 V / 




^IVie V / CJ 


^IVie V / C J 


^IVie V / CJ 


^IVie V / CJ 






pBe 


1 9Q _L_ /I 

Izo ± 4 


1 of- i r 

loo ± o 


lie i i o 

lib ± IZ 


loU ± lo 








13b ± 3 


146 ± o 


149 ± 25 


i or* i i n 

130 ± 19 


1 Q 9 ") 

±.y - z.o 


t^t" 1 — Trip ^ 0.27 


pru 


lol it o 


1 /|« _|_ Q 

141:0 it y 


1 9fi -1- 1 R 
IZO It ID 


1 9fi 4- 9^ 
1ZD It ZO 






DO 


1 /|Q _|_ A 


oon 4- osj 

ZZU T ZO 


017 4- oo 
Zl / ± Zo 


1 _i_ /| o 
100 It 4Z 






orO 


1 4- ^ 

iyo it o 


97D 4-1/1 
Z/U x 14t 


9D7 -1-17 
ZU 1 1 11 


1 87 4- 9fi 
15/ It ZO 






pBe 


Loo ± 4 


loU ± o 


izb ± b 


14z it 






pS 


170 ± 30 


148 ± 4 


131 ± 10 


143 ± 9 


2.3 - 2.9 


— m p < 0.68 


pPb 


195 ± 5 


157 ± 5 


141 ± 8 


140 ± 11 






SS 


210 ± 4 


238 ± 9 


190 ± 10 


171 ± 15 






SPb 


256 ± 4 


246 ± 4 


205 ± 7 


202 ± 8 



TABLE IV. Inverse slopes (T) extracted from fits of the data to Equation [|. The errors are 
statistical. Systematic errors are shown in Table ||. Also shown are the inverse slopes extracted 
from RQMD, version 1.08, after correction for weak decay feed-down. For RQMD, there is no 
significant difference in the inverse slopes if rope formation is turned off. 
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y 


System 


Protons 


Antiprotons 




pBe 


0.293 ± 0.008 ± 0.026 


0.060 ± 0.003 ± 0.005 




pS 


0.377 ± 0.009 ± 0.034 


0.080 ± 0.008 ± 0.007 


1 Q 9 Q 

i.y - z.o 


r, PA 


U.4Z0 ± U.U1U ± U.U4o 


u.u/u ± U.UUD ± U.UUrf 




Q Q 
DO 


A 07 4- H 1 K 4- M /1£ 


n \~\ ° i n 4- n hqq 
U.41o ± U.Uob ± U.Uoo 




oru 


1 q £ _l n _Lin 

lo.b ± u.4 ± i.y 


n vm 4_ n n/iQ 4_ n mc 
U. rfyi ± U.U4o ± U.lUo 




pBe 


0.158 ± 0.008 ± 0.014 


0.053 ± 0.006 ± 0.005 




vS 


0.204 ± 0.013 ± 0.018 


0.057 ± 0.006 ± 0.005 


2.3 - 2.9 


pPb 


0.251 ± 0.010 ± 0.025 


0.076 ± 0.007 ± 0.008 




SS 


4.51 ±0.20 ±0.42 


0.505 ± 0.038 ± 0.047 




SPb 


12.0 ±0.3 ±1.6 


1.10 ±0.05 ±0.15 



TABLE V. dN/dy for protons and antiprotons with statistical and systematic errors. 
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y 


System 


I 

Protons 


lopes 

Antiprotons 


No 
Protons 


Ropes 

Antiprotons 


1.9 - 2.3 


pBe 
pS 

pPb 
SS 

SPb 


0.211 ±0.005 

n oTr \ i a ai n 

0.279 ± 0.010 
0.433 ±0.019 
5.64 ±0.56 
16.66 ± 0.73 


0.033 ± 0.002 
0.039 ± 0.004 
0.057 ±0.007 
0.712 ±0.122 
1.050 ±0.105 


6.19 ±0.94 
16.2 ± 1.2 


0.221 ±0.065 
0.355 ±0.163 


2.3 - 2.9 


pBe 
pS 

pPb 
SS 

SPb 


0.179 ±0.003 
0.221 ±0.005 
0.295 ± 0.009 
5.252 ±0.155 
12.71 ±0.13 


0.048 ± 0.002 
0.058 ± 0.003 
0.095 ± 0.005 
1.051 ± 0.073 
1.589 ± 0.049 


4.98 ± 0.24 
11.22 ±0.18 


0.366 ± 0.064 
0.462 ± 0.037 



TABLE VI. Predictions for dN/dy for protons and antiprotons from RQMD, version 1.08. The 
contribution from weak decay feed-down is included in these calculations. 
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Cherenkov Beam Counter 
Beam Veto Counter 



FIGURES 

TO Counter 

Silicon Multiplicity Counter 




Beam direction 



MWPC3 




MWPC4 



UCAL 



D: Dipole Magnet 



C: Threshold Gas Cherenkov Counter 
H: Scintillator Hodoscope 
MWPC: Multi-Wire-Propotional Chamber 



Q: Superconducting Quadrupole Magnet UCAL: Uranium-Cu-Scintillator Calorimeter 



FIG. 1. The NA44 experimental set-up. 
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2 3 4 

Rapidity 



FIG. 2. The p or p acceptance in y and pr- The (f> acceptance (not shown) decreases from 2ir 
at p T = to 0.012 at p T = 1.6 GeV/c. 
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p 




Mass 2 (GeV 2 ) 



FIG. 3. Mass-squared distribution from Hodoscope 3 after pions have been vetoed by the 
Cherenkovs. 
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• Direct p □ p +A+Z a Direct p £ p+A+Z 
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FIG. 4. Invariant cross-sections as a function of tut— m p for central SPb collisions from RQMD, 
with (open symbols) and without (solid symbols) feed-down from weak decays. 
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• SPb □ SS A pPb & pS * pBe 
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FIG. 5. Invariant cross-sections for protons as a function of tut - m. 
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• SPb □ SS ▲ pPb -k pS O pBe 
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FIG. 6. Invariant cross-sections for antiprotons as a function of tht - rn p . 
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NA44Data □ RQMD with Ropes 
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^ * ff 



Systematic Error I 
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FIG. 7. Inverse slopes of the transverse mass distributions for each system for data and RQMD. 
Statistical and systematic errors for the data are added in quadrature. The global systematic errors 
common to all systems (Table ||) are shown by bars near the bottom right hand corner of each 
plot 
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NA44 Data □ RQMD with Ropes A RQMD no Ropes 
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FIG. 8. Rapidity densities (dN/dy) for each system for data and RQMD. Statistical and sys- 
tematic errors for the data are added in quadrature. 
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NA44Data □ RQMD with Ropes A RQMD no Ropes 
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FIG. 9. The integrated p to p ratio for data and RQMD. Systematic and statistical errors have 



been added in quadrature. The value of the integrated ratio for pp collisions [15] is 0.33 ± 0.13 



32 



20 
17.5 

15 
12.5 

10 

7.5 
5 

^2.5 

Z o 

Q 2 
1.75 
1.5 
1.25 
1 

0.75 
0.5 
0.25 




~T 



NA44Data □ RQMD Ropes A RQMD No Ropes 







Systematic Error I 



y=1.9-2.3 



4 



□ 



□ 



Systematic Error I 



y=2.3-2.9 



□ 



□ 



A 



A 



11-6% 



6-0% 11-6% 

Centrality 



6-0% 



FIG. 10. Proton and antiproton yields as a function of centrality for SPb collisions. Statistical 
errors are shown for each point. The vertical bar shows the systematic error of 13.6% common 
to all of the data, see Table ||. Also shown are the RQMD predictions for the same centrality 
selection. 
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FIG. 11. Invariant cross-sections as a function of mx — m p for central SPb collisions from data 
and RQMD. The RQMD distributions include the effect of feed-down from weak decays. 
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• NA44 □ NA35 
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FIG. 12. Comparison of NA44 and NA35, |iqjl7| ], rapidity densities for p — p and p. The pA 
data were taken at 200 GeV/c for NA35 and 450 GeV/c for NA44. The NA44 SS and SPb data 
have been selected to have approximately the same centrality selection as the NA35 SS and SAu 



data, and all NA44 data have been corrected for feed-down using the factors in Table [II. The 
NA35 p — p data are from the rapidity intervals 2.0-2.5 and 2.3-3.0 for both SS and SAu; for pA 
the NA35 rapidity ranges are 1.8-2.2 and the average of two bins covering rapidity 2.2-3.0. The 
NA35 p data are in the rapidity range 3-4. 
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